Origin of the anomalous heat current in coUisional granular fluids 
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We present a heuristic explanation for the anomalous (density-gradient-dependent) heat current 
in coUisional granular fluids. Inelastic grain collisions lead to highly non-equilibrium states which are 
characterized by large spatial gradients and/or temporal variations in the granular temperature, ft 
is argued that the heat current in such non-equilibrium states is driven by the temperature gradient 
averaged over the typical collision time and/or mean free path. Due to the density dependence of 
the inelastic energy loss, the nonlocal averaging of the temperature gradient leads to an effective 
dependence of the heat current upon the density gradient. 

PACS numbers: 45.70.Mg, 81.05.Rm, 44.10.-|-i 
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When energy is supplied to a system of inelastic grains, 
a particularly simple coUisional state can occur in which 
only binary collisions between the grains are important. 
The coUisional state occurs for hard grains at any den- 
sity below jamming, and for more realistic soft grains at 
densities sufficiently low that enduring and multiparti- 
cle contacts are rare 0, Q • In recent decades coUisional 
granular fluids have served as a paradigm for simulations 
and for the application of statistical ideas to granular 
systems 0, One perplexing phenomenon long pre- 
dicted for coUisional granular fluids is a violation of 
Fourier's law, which states that the heat current is pro- 
portional to the temperature gradient and always flows 
from hot to cold. In contrast, the heat current predicted 
by kinetic theory i i8j for a granular Uuid is 



-kVT — /iVn 



(1) 



where T = m{v'^)/d is the granular temperature in 
d space dimensions for particles of mass m and mean 
squared velocity (w^), and n is the grain number density. 
Here k is the granular thermal conductivity and ^ is a 
new transport coefficient that is nonzero for inelastic sys- 
tems. Due to the last term in Eq. ^ the anomalous heat 
current, it is possible when Vn ^ for the total heat cur- 
rent to Uow from cold to hot. Although counterintuitive 
the anomalous heat current does not violate the second 
law of thermodynamics, as fiuidized granular media are 
highly non-equilibrium systems. 

In the presence of gravity the anomalous heat current 
is manifested by a nonzero heat current at the height 
at which the temperature gradient is zero |^. More 
strikingly, for a system with a free upper surface (con- 
fined by gravity and excited from below), the anomalous 
heat current can lead to a temperature inversion: de- 
spite the strictly upward heat current, the temperature 
first falls as a function of height and then rises again in 
the vicinity of the free upper surface T^. In a recent 
experiment using NMR to probe vibroUuidized mustard 
seeds, Huan et al. ^ observed such a temperature inver- 
sion and quantitatively fit the heat current (including the 



anomalous term) to theoretical predictions. Fig. ^ shows 
example data from this study along with the hydrody- 
namic fit. A temperature inversion may also have been 
observed in much earlier experiments on vibroUuidized 
systems ^1 j but the inversion was not associated at 
the time with the anomalous heat current. 

Although a density-gradient-dependent heat current 
has been demonstrated in numerical simulations of gran- 
ular fluids 0,01, and ^ has been calculated using sophis- 
ticated statistical-physics methods 0, IS B (B] j a simple 
heuristic explanation of this phenomenon has not been 
presented to date. In this paper we present such an ex- 
planation, tracing the anomalous heat current directly to 
the highly non-equilibrium nature of coUisional granular 
fluids. Due to inelastic energy losses, such fluids necessar- 
ily have large time derivatives and/or spatial gradients of 
the granular temperature. As the temperature may vary 
significantly over a mean-free-path distance or grain col- 
lision time, the inherent nonlocality of the grain collision 
process must be taken into account. By considering the 
consequences of this nolocality, we show that the den- 
sity dependence of the inelastic energy loss leads to the 
anomalous heat current term. We argue that the to- 
tal heat current is fundamentally driven by temperature 
gradients, not density gradients, and that the apparent 
dependence of the heat current upon density gradients 
arises from re-expressing nonlocal quantities in terms of 
local ones. 

The existing technical calculations 0, S of the 
anomalous transport coefficient /i are valid for wide 
ranges of the dimensionless density na^ and the inelas- 
ticity parameter 1 — a (here a is the grain diameter and a 
is the normal restitution coefBcient). To validate our ex- 
planation for the anomalous heat current, we show that 
it predicts a value for fi that agrees up to factors of or- 
der unity with these more rigorous calculations in the 
dilute, nearly elastic limit na^ <^ 1, (1 — a) ^ I. The 
aim of the present paper is not to reproduce precisely the 
earlier results, but rather to provide a simple, physically- 
motivated explanation for the last term in Eq. ^ 
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FIG. 1: Measured granular density (top) and temperature 
(bottom) as functions of height for a dilute granular fluid 
that is vibrofluidized from below, from Ref. pH . The sym- 
bols show NMR measurements for mustard seeds vibrated at 
15g and 50 Hz, while the curves show a fit to the hydro- 
dynamic theory of Ref. Dimensionless units z" — z/a, 
n* = na^, T* = T/mga are used. The fitted restitution coef- 
ficient was a = 0.87. The sample was confined by gravity to 
the bottom region of the sample chamber with a free upper 
surface, so there was no vibrational energy input from above. 
Nevertheless, the experimental data show a clear minimum 
in the granular temperature T* as a function of height z* . 
At heights above the temperature minimum, the heat current 
(always upward) is from cold to hot. 



We begin with the commonly-used relaxation-time ap- 
proximation for the evolution of the heat current: 
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(2) 



where Vq = 



[T/mY/'^ is the relaxation rate for 
heat currents, roughly equal to the grain collision fre- 
quency. Here and below, hq,bu ■ ■ ■ are dimensionless con- 
stants of order unity that depend only upon the dimen- 
sion of space d; values are given in Table U Note that 
(T/to)^/^ is a typical particle speed and a'^^^ is roughly 
the collision cross section. 

The first term on the RHS of Eq. |21 gives the buildup 
of the heat current due to the free streaming of par- 
ticles from nearby regions with differing temperatures, 
while the second term expresses the relaxation of q to- 
wards zero due to collisions, which randomize the par- 



TABLE I: Dimensionless prefactors used in kinetic theory 
expressions in this paper, for dimensions d = 2 and 3. The bq 
values are calculated for a Maxwell-Boltzmann velocity dis- 
tribution, expected to be a good approximation for weakly in- 
elastic systems fl3|. The hi, and 6^ values are computed from 
the results tablulatcd in Appendix A of Ref. . Finally, the 
b\ values are those required for the anomalous transport co- 
efficient /i to be the same for cooling states and steady states 
according to the approximate calculation presented in this 
paper. 
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Cooling rate 




80F/3 
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Mean free path 
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tide velocities. In an elastic system {a — I) perturbed 
by weak temperature and density gradients, the temper- 
ature varies slowly in time due to energy conservation. 
Hence a steady state is reached with dq/dt « 0, giv- 
ing Fourier's law: q = — kqVT with kq = bqiiT /mvq ^ 
(r/TO)i/Va'*-i. 

Implicit in this derivation of Fourier's law are two re- 
lated "near-equilibrium" assumptions: (i) dynamics af- 
fecting n and T occur very slowly relative to the heat 
current relaxation rate Vq^ and (ii) spatial variations of n 
and T occur on characteristic length scales much greater 
than the mean free path \q ~ \/na'^~^. In an inelas- 
tic system (a < 1) the near-equilibrium assumptions are 
not generally valid because the system continually loses 
energy, which causes temperature changes at rates com- 
parable to Vq. If energy is supplied at the boundary of the 
system to maintain a steady state, this leads to density 
and temperature variations on length scales comparable 
to Xq. To account for these non-equilibrium effects the 
derivation of Fourier's law must be generalized to be non- 
local in space and time. We show here how this leads to 
an effective dependence of the heat current upon Vn. 

As a first example, we consider a. freely- cooling granular 
fluid (i.e., no energy input) with small gradients in T and 
n. Locally, one has dT/dt ~ ~C,T with the cooling rate 
C given to leading order in the inelasticity (1 — a) by: 
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(3) 



For typical values of the restitution coefficient a w 0.9, 
the ratio of the cooling rate to the heat current relaxation 
rate, C,/^q = (^c/^!^)(l ~ ^ small but not infinites- 

imal quantity (w 0.1 in 3D and 0.3 in 2D). Therefore 
T changes non-negligibly during a typical grain collision 
time l/fc'g, which violates the first near-equilibrium as- 
sumption above. 

We use the following simple argument to estimate the 
heat current q(i) in this situation. According to Eq. [21 
q(t) continuously relaxes towards the Fourier-law value 
— kqVT, which however is time-dependent. Thus, q(t) 
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lags behind the Fourier-law value by approximately one 
relaxation time ^/vq- 



q(<) « -KaVT + {l/vq)d{HiQVT)/dt. 



(4) 



Equation ^ expresses the idea that the gradient that 
drives the heat current is nonlocal in time: The heat 
current at time t reflects the weighted average of —kqVT 
over earlier times, with an typical lag of one collision time 
The second term in Eq. ^is readily calculated, 



diKoVT)/dt = ^VT 
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where C cx nT^/^ and Kg oc T^/"^ have been used to sim- 
plify this equation. Thus: 
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From Eq. |21we can read off the anomalous transport co- 
efficient /X = kqCT /vqU, proportional to the semi-small 
parameter C,/vq. The ordinary thermal conductivity also 
acquires a correction proportional to the same parame- 
ter. Both of these results are close to the results of more 
sophisticated calculations [E IE 0) for example the 
value of /i is 5 /4 times the result calculated in Ref . . 

The origin of the heat current's effective density- 
gradient-dependence is the dependence of the cooling rate 
on density {dC,/dn ^ in Eq.|SJ). For example, if at time 
t = there is a density gradient but no temperature 
gradient, then at earlier times the denser regions must 
have been hotter since they are cooling faster. There- 
fore, averaging kqVT over times t < Q yields a non-zero 
heat current in the absence of a temperature gradient at 
t = 0. 

As a second example we consider a boundary- 
driven steady-state — e.g., a vertically- vibrated granular 
medium such as that represented in the data of Fig. ^ 
and described in Ref. [ij. Due to inelastic energy losses 
and gravity, a boundary-driven steady-state necessarily 
has strong spatial variations of n and T, in violation of 
the second near-equilibrium assumption listed above. 

A rigorous treatment of transport in the presence of 
strong gradients is complex, but the folowing simple ar- 
gument should give a reasonable approximation. The 
typical distance over which a particle drifts before hav- 
ing its velocity randomized by a collision is the mean free 
path \q — h\{T jraf-l'^ jvq ^ l/n(T'^~^, where bx is an- 
other dimensionless constant of order one. Therefore we 
should replace VT in Eq.|21by (VT)a, which denotes the 
spatial average of VT over the mean free path Xq. Thus, 
in the steady state the temperature gradient is effectively 
nonlocal in space, just as it is effectively nonlocal in time 
for the cooling state. 




FIG. 2: Schematic illustration of the mechanism for anoma- 
lous heat currents in a steady-state situation similar to that 
shown in Fig. At the height z of the temperature mini- 
mum, the gradient dT/dz has non-zero curvature due to the 
density gradient and the density dependence of the inelastic 
power loss. As a result of the curvature, the average of the 
temperature gradient over the mean free path Aq is non-zero 
(filled dot on graph). This non-zero effective gradient in turn 
drives a heat current even though dT/dz is zero at this height. 
(Although it might appear from this diagram that an anoma- 
lous heat current can only occur over a small range of z near 
the minimum, this is not the case due to the increase of Xq 
with z.) 



In the steady state {VT)\ differs from VT by an 
amount proportional to the curvature of VT times A^, 
as shown schematically in Fig. |2| To simplify the calcu- 
lation, we assume that T and n depend only upon the z 
coordinate (typically in the direction opposite gravity). 
Then 
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dT/dz+{Xl/&)dT^/dz^. (7) 



To lowest order in C,/vq we have dT/dz = —qz/i^f), where 
qz is the z-component of the heat current q. Hence 



dT'/dz' « --K-'d'qz/dz 



(8) 



This omits corrections proportional to derivatives of kq, 
which considerably complicate the calculation but do not 
affect the leading-order estimate for /z. 

To complete the argument we note that in a steady 
state the energy loss rate per unit volume due to inelastic 
collisions, Pc, is balanced by the negative divergence of 
the heat current: — — V ■ q = —dq^/dz. The energy 
loss rate is proportional to the cooling rate C,: Pr = cC/T = 
dnC,T/2 where c = dn/2 is the specific heat Using 
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Eqs. |7|and|Sl we have 

{dT/dz)x - dT/dz 



-{Xy&>io)d\Jdz^ 
= {Xld/12Ko)d{nCT)/dz 
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where C oc nT^/'^ has again been used. Thus we find for 
the steady-state heat current: 



qz = -Ko{dT/dz)x 
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giving /i = b\dKQ<^T l&hqVqTi for the anomalous heat 
transport coefficient in a boundary-driven steady-state. 
This value for /i in the steady state agrees with the 
cooling-state estimate of derived above for b\ = 
(66q/(i)^/^, which gives reasonable values for the dimen- 
sionless constant bx (Table QJ. 

For the boundary-driven steady-state, the origin of the 
anomalous heat current is again the density dependence 
of the coUisional energy loss (Pc = cC,T ~ v?), which 
forces the heat current and hence the temperature gra- 
dient to have a large, non-zero second spatial derivative 
(curvature). The curvature of the temperature gradient 
over distances comparable to the mean free path in turn 
implies that the effective (nonlocal) temperature gradi- 
ent (VT)a is different from the local gradient, leading to 
an extra term in the heat current that is proportional to 
the density gradient (Fig.|21). 



In summary, although the anomalous heat current is 
nominally proportional to the density gradient, an analy- 
sis of the underlying physical effects shows that the driv- 
ing term is actually the temperature gradient. Due to 
the large time and/or space derivatives of the tempera- 
ture that necessarily occur in dilute granular fluids, the 
driving temperature gradient must be averaged over a 
typical collision time l/vq and mean free path Ag. Thus 
the effective gradient is nonlocal in time and space; when 
expressed in terms of local gradients an apparent density- 
driven heat current appears. 



Finally, we emphasize that it is the strong space or time 
derivative of the temperature, rather than the inelasticity 
per se, that is the source of the anomalous heat current. 
Therefore, similar heat currents should in principle be 
observable in elastic systems if suitably strong tempera- 
ture gradients could be imposed. From this point of view, 
the role of inelasticity in a granular fluid is primarily to 
create a highly non-equilibrium state, and some novel 
physical properties of granular fluids (like the anomalous 
heat current) are due fundamentally to the large devi- 
ation from equilibrium rather than the inelastic energy 
loss. 
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